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I.  INTRODUCTION 


The  possible  applications  of  railguns  for  military  uses  have  been  given 
considerable  attention  in  recent  years.  Of  the  many  classes  of  electric  guns, 
the  simple  linear  railgun  has  received  the  most  study.  Within  this  class,  by 
far  the  most  successful  has  been  the  plasma  armature  type  exemplified  by  the 
CHECMATE  railgun.^  In  this  device,  a  thin  aluminum  foil  on  the  rear  face  of 
the  projectile  is  subjected  to  a  large  current  and  rapidly  vaporizes.  That 
current  Is  Initiated  in  the  circuit,  which  includes  a  substantial  inductor, 
by  the  discharge  of  a  capacitor.  At  the  Instant  of  complete  discharge,  the 
capacitor  is  then  shunted  out  of  the  circuit  by  a  low  resistance  "crowbar". 
Driven  by  the  energy  stored  in  the  field  around  the  inductor,  the  current  con¬ 
tinues  to  flow  through  the  vapor  left  by  the  explosion  of  the  foil,  heating 
and  ionizing  it  to  produce  an  arc  or  plasma.  The  driving  force  applied  to  the 
projectile  arises  from  the  pressure  accrued  in  the  gaseous  plasma  as  the  dif¬ 
fuse  current  of  density  J  flowing  through  it  from  rail  to  rail  reacts  with  the 
accompanying,  intense  magnetic  induction  field,  B.  The  term  plasma  (or,  altej*- 
natively,  arc)  implies  a  hot,  ionized,  highly  luminescent,  gaseovis  conduction 
path  for  the  railgun  c\rrrent;  the  term  armature  relates  to  the  motive  part  of 
an  electric  motor.  Whfle  the  exact  environment  between  the  rails  behind  the 
projectile  is  not  known  in  great  detail,  several  theoretical  descriptions^"^ 
of  its  state  have  been  offered  and  numerous  experimental  studies  have  been 
performed  on  various  small  railguns. 

Most  of  the  projected  uses  of  railguns  do  not  Involve  small  devices,  but 
rather  railguns  whose  bore  dimensions  are  several  centimeters  and  whose  lengths 
are  several  meters.  The  CHECMATE  railgun  has  a  square  bore  5  cm  by  5  cm  and  a 
length  of  five  meters.  It  has  been  used  to  fire  projectiles  having  masses 
between  100  and  350  grams  at  velocities  between  one  and  three  kilometers  per 
second.  At  the  time  of  this  writing,  it  has  been  fired  nearly  100  times.  The 
successful  operation  of  the  CHECMATE  railgun  represents  a  significant  advance 
toward  practical  applications  of  railgun  technology.  A  schematic  illustration 
of  the  railgun  and  its  power  supply  are  shown  in  Fig,  1, 

In  this  report,  we  describe  a  set  of  diagnostic  meastirements  made  on  the 
CHECMTE  railgun  in  an  attempt  to  quantify  some  of  the  properties  of  the  plasma 
arc  armature.  In  addition,  these  measurements  are  used  in  conjimction  with  a 


^Holland,  M.M. ,  Wilkinson,  G.M. ,  Krickhuhn,  A.P,  and  Dethlefson,  R. ,  "Six 
Megajoule  Rail  Gun  Test  Facility,"  IEEE  Trans,  Magn.  MAG-22,  1521  (1986). 

^Powell,  John  D.  and  Batteh,  Jad  H.,  "Plasma  lynamlcs  of  an  Arc-Driven  Elec¬ 
tromagnetic  Projectile  Accelerator,"  J,  Appl,  Phys,  52.  2717  (1981).  See 
also,  "Plasma  Dynamics  of  the  Arc-Driven  Rail  Gun,"  Ballistic  Research  Labora¬ 
tory  Report  No.  ARBRL-TR-0226T,  September  1980. 

3powell,  John  D.  and  Batteh,  Jad  H,,  "Two-Dimensional  Plasma  Model  for  the 
Arc-Driven  Railgun,"  J.  Appl.  Phys.  5^,  22U2  (1983).  See  also  Powell,  John 
D, ,  "Two  Dimensional  Model  for  Arc  Dynamics  in  the  Railgun,"  Ballistic  Re¬ 
search  Laboratory  Report  No.  ARBRL-TR-02i*23,  October  1982 

^Sloan,  M.L,  "Physics  of  Rail  Gun  Plasma  Armatures,"  IEEE  Trans.  Magn.  MAG-22 , 
17h7  (1986). 

^Thio,  Y.C.  and  Frost,  L.S.,  "Non-Ideal  Plasma  Behavior  of  Railgun  Arcs,"  IEEE 
Trans.  Magn.  MAG-22.  1T57  (1986). 
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Figure  1.  Schematic  Illustration  of  CHECMATE  Railgun  and  Power  Su’ 


one-dinensional,  steady-state  model  to  infer  certain  properties  of  the  arc 
which  cannot  he  measured  directly.  The  report  is  organized  as  follows.  In 
Sec.  II  we  describe  the  diagnostic  probes  constructed  especially  for  these 
experiments,  indicate  the  specific  measiorements  made,  and  analyze  the  data  to 
infer  the  current  density  and  arc  length.  In  Sec.  Ill,  we  indicate  some  re¬ 
visions  undertaken  on  a  model  previously  developed  by  Powell  and  Batteh,^ 
and  then  use  the  revised  model  to  calculate  properties  of  the  plasma  that 
cannot  be  measured  directly.  In  Sec.  IV  are  contained  our  conclusions  and 
recommendations  for  further  study. 


II.  EXPERIMENTAL  RESULTS 

1.  Diagnostic  Probes 

The  measurement  of  electrical  qmntities  on  a  railgun  requires  particular 
care  to  minimize  the  interactions  of  its  characteristically  intense  electric 
and  magnetic  fields.  Probes  must  be  provided  with  the  maximum  practicable 
insulation  possible  and  multiple  grounds,  establishing  loops  in  which  currents 
may  be  magnetically  induced,  must  be  diligently  avoided.  All  of  the  diagnos¬ 
tic  techniques  employed  have  been  previously  used  and  refined  on  the  small 
BRL  railgun.®  The  voltages  produced  by  the  probes  were  digitally  sampled  at 
a  preset  rate  using  Ilicolet  oscilloscopes.  These  histories  were  then  stored 
on  floppy  discs  after  each  firing.  The  discs  were  later  read  into  the  BRL- 
LFD  VAX  8600  for  analysis.  Each  type  of  probe  and  its  response  characteris¬ 
tics  are  briefly  described  in  the  sections  following.  The  only  additional 
data  used  in  this  analysis  were  railgun  cxirrent-time  histories  and  the  shot- 
start  locations.  To  obtain  the  gun's  current-time  histories,  signals’^  were 
recorded  from  calibrated  Rogowski  loops  encircling  each  of  the  five  power- 
supply  module-output  leads.  These  signals  were  integrated  and  summed  to  give 
the  total  overrent  delivered  to  the  railgun. 

The  conditions  and  performance  parameters  of  each  of  the  four  firings  are 
s\immarized  in  Table  I.  For  reference  we  shall  utilize  the  shot  n\imber  from 
the  Maxwell  firing  log. 


TABLE  I.  Description  of  Shot  Parameters 


Maxwell 

Shotf^ 

Projectile 

Mass 

Peak 

Current 

Velocity 

Injection 

Velocity 

Capacitor  Cl 
Voltage 

9h 

120  g 

1.210  MA 

2.37  km/s 

680  m/s 

35  kV 

95 

233  g 

1.2lh  MA 

1.86  km/s 

500  m/s 

35  kV 

96 

120  g 

1.215  MA 

2.36  km/s 

680  m/s 

35  kV 

97 

233  g 

l.hOk  MA 

2.105  km/s 

500  m/s 

39  kV 

^Jamison,  K.A.  and  Burden,  H.S.,  "A  Laboratory  Arc  Driven  Rail  G\in,"  Ballis¬ 
tic  Research  Laboratory  Report  No.  ARBRL-TR-02502,  June  1983. 

"^Wilkinson,  G.M.  (unpublished  data). 


a.  Voltage  Probes 


Since  very  different  voltages  exist  at  the  two  ends  of  the  railgun, 
the  breech  end  muzzle  ends  of  the  rails  were  each  shunted  by  a  resistive  chain 
that  was  j.inked  through  a  terminated  Pearson  current  transformer  (see  Fig.  2); 
the  voltage  output  of  the  Pearson  was  thus  proportional  to  the  source  voltage. 
Series  diodes  Isolated  the  probe  circuits  from  the  discharge  of  CHECMATE's 
negatively  charged  sense  capacitor  on  passage  of  the  projectile. 

b.  B-Dot  Loops 

In  these  measiirements ,  the  time  rate  of  change  of  the  armature  induc¬ 
tion  field  is  sensed  in  a  coll  oriented  to  be  insensitive  to  the  fields  pjro- 
duced  by  the  ciirrents  in  the  rails.  Essentially,  the  coil  axis  must  be  para¬ 
llel  to  the  rail  currents  and  pjerpandicular  to  the  curmature  current.  When  an 
armature  carrying  megamjare  ciorrents,  with  accompanying  10  to  25  Tesla  induc¬ 
tion  fields,  passes  a  point  at  speeds  of  the  order  of  1  km/s,  the  rising  and 
falling  field  intensity  can  develop  a  significant  induced  voltage  in  even  a 
small  conducting  loop.  The  relative  orientations  and  typical  output  for  such 
a  coil  are  shown  schematically  in  Fig.  3  together  with  the  governing  equations 
for  the  induced  signals.  Table  II  denotes  the  locations  and  dimensional  de¬ 
tails  for  each  of  the  stations  mounted  in  the  CHECMATE  barrel  for  these  exp)erl- 
ments.  With  five-turn  coils  less  than  one  centimeter  in  diameter,  we  have 
measiured  inductive  signatures  up  to  seven  volts  in  amplitude  \ising  the  probe 
assembly  shown  in  Fig.  1*.  This  assembly  was  designed  particularly  for  com¬ 
pactness  so  that  the  hole  in  the  sidewall  of  the  CHECIiATE  barrel  needed  to  ac¬ 
commodate  it  would  cause  the  least  possible  weakening  of  the  structure.  The 
modification  of  the  CHECMATE  sidewall  to  accept  the  probe  is  shown  in  Fig.  5. 
The  view  of  the  probe  shown  in  Fig.  It  is  rotated  90°  with  respject  to  the  per¬ 
spective  in  Fig.  5. 


TABLE  II. 


Probe  Stations  and  Dimensions 


Distance  from  Arc  Initia¬ 
tion;  Shots  9^  and  96 
Distance  from  Arc  Initia¬ 
tion;  Shots  95  and  97 
Number  of  Turns 

Center  of  Bore  to  Loop 
Center 

Loop  Diameter 


Station  1 
58.it  cm 

68.6  cm 
5 

5.30  cm 
0.876  cm 


Station  2 

83.8  cm 

9it.O  cm 
5 

5.36  cm 
0.876  cm 


Station  3 
109.2  cm 

119«i*  cm 
5 

5.33  cm 
0.876  cm 


Station  it 
13it.6  cm 

litit,8  cm 
5 

5.it0  cm 
0.876  cm 


•NOTE:  The  exact  location  of  arc  initiation  varies  with  the  injection  veloc¬ 
ity  which  is  different  for  each  projectile  mass. 


c.  Flber-Ontic  Light  Probes 


The  plasma  armature  is,  to  understate  the  case,  highly  luminous. 
Our  experience  in  routing  light  from  the  bore  to  a  reverse  biased  PIN  diode 
has  indicated  that,  even  without  drilling  completely  through  the  insulating 
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Figure  2.  Schematic  Circuit  for  Breech  and  Muzzle  Voltage  Probes. 
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Figure  4.  Induction  Field  and  Light  Emission  Probe  Assembly. 
A  cross  section  is  shown;. 


CO 

C/)  _) 

2  a:  _j 

QL  ZD  u  _]  a: 

IjJ  -J  CD  O'  CD 

Q. _]  o  q:  u 


2  U 

cr  u 
u 

U  CkC 
CD 

o 

CkC 

0. 


8 


Figure  5.  Cross  Section  View  of  Probe  Mounting  Location  in  CHECMATE  Barrel  Structure 


sidewalls  of  the  railgun,  a  satisfactory  signal  will  he  generated  as  the 
plasma  armature  passes.  Figures  ^  and  5  show  the  method  of  positioning  the 
fiber  optic.  The  thickness  of  insulator  separating  the  bottom  of  the  fiber 
optic  cavity  and  the  bore  surface  was  0.6U  cm;  the  choice  of  this  thickness 
represents  some  conservatism  over  that  originally  specified. 

The  fiber  optics  were  routed  to  a  well-shielded  box,  physically  and  elec¬ 
trically  separated  from  the  gun  barrel,  said  which  enclosed  the  PIN  diodes 
used  to  sense  the  light.  The  remainder  of  the  circuit  is  shown  in  Fig.  6. 

The  battery  was  connected  with  polarity  such  that  the  diode  blocked  current 
flow.  As  the  diode  is  illuminated,  however,  liberated  electrons  constitute  a 
current  proportional  to  the  incident  light  flux  and  the  voltage  drop  across 
the  diode.  This  current  produces  a  voltage  drop  in  the  load  resistor  nearly 
proportional  to  the  li^t  intensity.  Because  the  a]!^>lltude  of  this  voltage 
varies  over  a  range  limited  by  the  supply  voltage,  the  original  9”Volt  battery 
was  changed  to  67  volts  to  con5)ensate  for  the  increased  cavity  wall  thickness 
discussed  above. 

2.  Measured  Quantities 


The  probes  described  in  the  previous  section  were  utilized  to  meas\ire  the 
following  quantities: 

o  Rail  (armature)  current 

o  Breech  voltage 

o  Muzzle  voltage 

o  Time  derivatives  of  the  armature  induction  field 

°  Light  emission  to  the  insulating  portion  of  the  barrel 
o  Armature  location  versus  time. 

a.  Ciirrent 

Figure  7  is  a  representation  of  the  current-versus-time  plots  measured 
by  Integrating  and  slimming  the  Maxwell  Rogowskl  coil  signals  from  each  of  the 
five  capacitor  banks.  The  waveform  is  typical  of  a  crowbarred  inductor  cir¬ 
cuit  which  is  characterized  by  a  rapid  rise  to  peak  current  followed  by  a 
nearly  exponential  decay.  Current  waveforms  from  Shots  9^  and  96  are  nearly 
identical  and  are  indistinguishable  in  Fig.  7.  Table  I  above  lists  for  each 
of  the  four  shots  a  value  of  peeik  current.  In  practice,  for  the  times  when 
the  armature  is  near  the  probe  array,  the  cuirrent  may  be  accurately  modeled  by 


I(t)  =  1(0)  exp  (-t/t) 


(1) 


This  equation  allows  two  quantities,  l(0)  and  t,  to  be  used  in  the  computation 
of  the  Instantaneous  current  at  times  when  the  armature  is  close  to  the  probe 
array.  Table  III,  in  Sec.  e  below,  gives  values  for  both  the  scaling  con¬ 
stant  and  the  decay  constant  which  describe  the  railgun  current  for  all  four 
shots . 
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Figure  6.  Schematic  Circuit  of  PIN  Diode  Light  Detectors 


b.  Breech  and  Muzzle  Voltage 

The  voltage  across  the  extreme  tip  of  the  arc  drives  the  current 
sensed  in  the  muzzle  voltage  circuit.  This  circuit  contains  only  portions 
of  the  rails  in  which  no  propulsive  current  is  flowing  so  self-inductive 
and  resistive  effects  of  these  large  currents  are  not  superimposed  on  the 
measurement.  (Mutual  inductance  between  the  propulsive  and  the  measviring 
circuits  may  also  be  shown  to  be  very  small.)  The  voltage  seen  is,  then, 
the  resistive  voltage  across  the  leading  edge  of  the  arc.  The  breech 
measurement  path,  on  the  other  hand,  includes  both  ciarrent  cairyliig  rails 
and  the  trailing  edge  of  the  eirmature.  This  path  surrounds  a  growing  ' 
volume  of  flxix  driving  the  armature.  The  resulting  inductive  voltages 
are,  thus,  added  to  the  resistive  voltage  developed  across  the  rear  por¬ 
tion  of  the  armature.  Figure  8  shows  breech  and  muzzle  voltage  traces  for 
Shots  96  and  97.  The  rapid  fall  of  both  voltage  traces  dxnring  the  early 
portion  of  the  shot,  0  to  300  microseconds,  signals  the  formation  of  the 
armature  plasma.  After  500  microseconds  the  voltage  drop  becomes  relatively 
constant  at  or  near  300  V.  The  muzzle  voltage  of  only  300  V  is  lower  than 
anticipated  for  a  5  cm  bore  railgun.  Of  particular  interest  are  the  muzzle 
voltage  jumps  at  2.U  ms  and  2.9  ms  for  Shots  96  and  9T,  respectively.  These 
jumps  are  characteristic  of  the  armature  exit  from  the  railgun  when  there 
is  still  sizable  field  energy  in  the  inductor  and  between  the  rails.  We  also 
note  the  lack  of  such  jumps  in  the  breech  voltage  suggesting  that  the  armature 
was  quite  long  near  the  end  of  the  shot. 

c.  Armature  Induction  Field 


The  armature  induction  field,  B,  is  not  directly  sensed;  instead,  its 
time  derivative,  "B-dot,"  is  measured  at  four  locations.  This  measurement 
constitutes  the  standard  against  which  a  modeled  cotinterpart  is  compared  as  a 
test  of  the  success  in  the  modeling  procediure  developed  to  unfold  the  armature 
cvirrent  density.  The  characteristic  form  of  a  B-dot  trace  is  a  relatively 
sharp  negative  peak  followed  a  wider,  lower  amplitude  positive  pesdc..  The 
polarity  of  these  peaks  is  tied  to  the  chosen  orientation  of  the  B-dot  coil. 

As  the  distributed  current  in  the  arc  approaches  the  array  of  coils,  each  of 
its  elements  is  contributing  to  each  respective  fall  in  voltages  sensed.  As 
the  midpoint  of  the  distribution  passes  a  given  coil,  the  effect  of  the  por¬ 
tion  of  the  current  which  has  already  passed  and  is  receding  matches  the  eff¬ 
ect  of  the  portion  still  to  come  and  the  response  passes  throu^  zero.  Finally, 
as  the  bulk  of  the  current  is  receding  down  bore,  the  polarity  reverses. 

Figures  9  throu^  12  show  the  B-dot  traces  for  each  of  the  four  stations 
for  Shots  9**  through  97.  The  indlvidtial  points  are  the  digitized ‘voltages  in¬ 
duced  in  the  loops  imbedded  in  the  barrel  and  smooth  lines  are  those  calculated 
by  the  fitting  procedure  to  be  discussed  in  the  next  section.  The  random 
point-to-point  variations  are  apparently  real  euid  the  disturbances  superimposed 
on  the  later  portions  of  the  traces  may  be  caused  by  small  arcs  behind  the 
original  one.  The  nearly  equal  amplitudes  from  the  different  probes  suggest 
that  the  increasing  length  of  the  arc  and  the  decay  of  the  total  cvirrent  which 
reduce  the  B-dot  signals  are  nearly  balanced  by  the  increase  in  velocity  which 
boosts  the  B-dot  signals. 
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d.  Light  Emission 

The  observed  intensity  of  emitted  light  as  a  function  of  time  at  each 
of  the  four  probe  stations  is  plotted  for  Shots  95  and  96  in  Figs.  13  and  14. 
The  differences  in  vertical  scale  maxima  are  an  indication  of  the  variation 
in  sensitivity  of  the  uncalibrated  PIN  diodes,  and  should  not  be  used  to  in¬ 
fer  relative  intensities  among  stations.  For  Shot  95,  the  initial  sharp  rise 
in  light  intensity  indicates  passage  of  the  armature-projectile  interface; 
for  Shot  96,  however,  luminosity  ahead  of  the  sharp  rise  indicates  leakage  of 
luminous  material  past  the  projectile  into  the  bore  ahead  of  it,  from  500  to 
700  microseconds  at  probe  station  2.  The  extended  leakage  path  and  greater 
stability  of  the  longer,  heavier  projectile  used  in  Shot  95  offer  an  explana¬ 
tion  of  this  observation.  On  Shot  94  our  digitizing  instrument  was  set  on 
too  sensitive  a  scale  and  the  signals  saturated  the  amplifiers.  On  Shot  97 
the  recording  equipment  failed  to  record  all  but  one  of  the  light  emission 
signals. 


e.  Armature  Location  Versus  Time 


Our  method  of  analysis  of  armature  current  density  from  B-dot  signa¬ 
tures  requires  knowledge  of  the  location  of  the  leading  edge  of  the  arc  ver¬ 
sus  time.  Normally,  we  would  use  the  sharp  rise  in  the  fiber  optic  data  to 
determine  the  position-versus-time  function;  however,  the  loss  of  some  fiber 
optic  data  and  the  uncertainty  introduced  by  the  precursor  luminosity  described 
above,  have  led  us  to  use  the  time  of  the  first  negative  peak  of  each  of  the 
four  B-dot  stations  to  determine  the  coefficients  of  the  cubic  equation 

Xp(t)  t’  +  A^  t^  +  A^  t  +  A^  .  (2) 

Here  Xp(t)  denotes  the  position  of  the  projectile-armature  interface.  From 
this  information  (see  Table  III),  the  acceleration  may  be  derived  and  plots 
of  acceleration  versus  time  are  shown  in  Figs.  15  and  16.  The  accuracy  of 
this  method  is  probably  only  on  the  order  of  10  percent,  but  it  is  good 
enough  to  show  that  the  ideal  accelerations,  predicted  from  the  effect  of  the 
calculated  armature  forces  on  the  projectile  mass,  are  too  high.  Bore  fric¬ 
tion  and  energy  required  to  drive  a  shock  ahead  of  the  projectile  could  ac¬ 
count  for  most  of  this  discrepancy. 


TABLE  III.  Shot 

Parmeters  1 

Used 

in  Data  Analysis 

Shot  94 

Shot  95 

Shot  96 

Shot  97 

1(0) 

1.30  Ma 

1.28  Ma 

1.30  Ma 

1.49  Ma 

T 

2.17  ms 

2.55  ms 

2.17  ms 

2.50  ms 

-1.1802  X  10® 

-4.9575  x 

io'7 

-7.254  X  10^ 

-3.4926  X 

10 

^2 

7.0107  X  10® 

3.9912  X 

10® 

6.5462  X  10® 

4.1452  X 

10 

A, 

944.27 

757.80 

938.40 

801.34 

-0.12780 


-0.14770 


-0.13405 


-0.11603 
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Figure  13.  Light  Intensity  Versus  Time  Traces  from  Four  Probes;  Shot  95 
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Figure  .  Light  Intensity  Versus  Time  Traces  from  Four  Probes;  Shot  96 
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Figure  16.  Ideal  and  Interpolated  Accelerations  for  CHECMATE  Railgun  Shots  95  and  97 
The  projectile  mass  was  233  g. 


3.  Data  Analysis 


a.  Current  Density  from  B-Dot  Signatures 


A  sovirce  of  further  understanding  of  the  plasma  armature  properties  may 
be  found  by  setting  the  time  derivative  of  the  Biot-Savart  integral  (Fig,  3) 
equal  to  instantaneous  values  of  the  field  derivative  sensed  by  a  B-dot  coll. 
Then,  an  unfolding  operation  may  reveal  the  kernel  of  the  integral  J{x,t)  which 
is  the  distribution  of  current  within  the  armature.  Because  there  is  no  reason 
to  believe  that  only  one  distribution  of  current  is  capable  of  producing  the 
observed  field  distribution,  not  enough  Information  is  available  to  allow  a 
unique  analytic  solution  to  this  problem.  In  ova*  previous  work,®  however,  we 
have  developed  a  computer  technique  for  modeling  a  moving,  current  distribu¬ 
tion  and  the  effect  of  its  magnetic  field  on  a  stationary  B-dot  coil.  We 
artificially  provide  added  information  by  using  a  modification  of  the  Gaussian 
cvirve  function,  truncated  at  the  plasma-projectile  interface  (see  Pig.  IT): 

J(x,t)  =  (Iq/Cj^)  exp  [-(x-X^)^/2W^] 

where  X  =  Xp-X(t).  Figure  IT  and  Eqs.  (3)  -  (5)  define  the  five  adjustable 
parameters,  P.. 


The  procedure  systematically  varies  the  parameters  until  the  calculated  coil 
response  closely  fits  the  B-dot  signatures  we  have  measured;  in  this  sense, 
the  distribution  is  "unfolded"  from  the  B-dot  signature. 


We  have  attempted  the  modeling  procedure  with  a  number  of  different  func¬ 
tions  but  were  most  successful  in  matching  the  set  of  B-dot  traces  when  using 
a  truncated  Gaussian  ciirve  function.  This  Gaussian  (or  statistical  normal) 
representation  of  the  current-density  distribution  was  chosen  for  several  rea¬ 
sons:  First,  it  contains  several  parameters  whose  variation  manipulates  it 
to  the  desired  shape.  Second,  the  shape  is  qualitatively  similar  to  standard 
theoretically  derived  profiles, ^  Third,  the  formality  for  normalizing  the  area 
is  well  known  and  easily  usable.  The  latter  requirement  is  in  accordance  with 
the  practical  necessity  that  the  area  (i.e.,  the  integral  of  the  current  den¬ 
sity  over  the  distribution)  should  be  normalized  to  the  value  of  the  instant¬ 
aneous,  total  current.  This  cvirrent  is  computed  using  the  values  of  t  and 
1(0)  (an  adjusted  initial  current  value)  and  the  given  value  of  time.  The 
value  for  1(0)  is  fitting  parameter  P5  smd  is  adjusted  to  improve  the  fit. 

The  value  for  t  is  fixed  at  the  experimentally  measured  quantity,  A  series 
of  calculations,  each  at  a  new  position  as  determined  from  the  digitizer  time 
step  and  Eq.  (2),  describing  this  iMving  distribution's  magnetic  field,  were 
made  using  a  Biot-Savart  voliune  integral.  The  time  derivative,  B-dot,  at  the 
position  of  a  given  coil  was  then  calculated.  With  each  trial  set  of  para¬ 
meters,  a  test  of  the  fit  between  the  calculated  B-dot  and  that  seen  by  the 
coil  was  made.  The  process  then  directed  appropriate  changes  in  the  peirame- 
ters  and  was  repeated  iteratively  until  no  further  improvement  in  goodness  of 
fit  could  be  obtained. 


®Jamison,  K.A.,  Marquez-Reines ,  M.  and  Burden,  H.S.,  "Measurements  of  the  Spa¬ 
tial  Distributions  of  Current  in  a  Rail  Gun  Arc  Armature,"  IEEE  Trans.  Magn. 
MAG-20.  i*03  (198^). 
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The  parameters  used  to  describe  the  current-density  distributions  which 
fit  the  data  best  are  listed  in  Table  IV.  Referring  to  Fig,  17,  we  see  that 
X(t),  the  distance  between  the  peak  of  the  Gaussian  and  the  front  tr\incation 
point  representing  the  back  of  the  projectile  is  given  by 


X(t)  =  P. 


/l{0)/4 

^TTtT^ 


=  V 


Pj,/^ 


(3) 


TABLE  IV.  Ciorrent  Density  Parameters  for  Acceptable  B-Dot  Pit 


Shot  9k 

Shot  95 

Shot  96 

Shot  97 

O.O692U 

0.07375 

O.O68OO 

O.O66U6 

^2 

0.03701 

0.01*986 

0.0391*5 

0.03699 

"3 

1.0717 

l.llOU 

1.2726 

1.1527 

-0.0902 

0.5518 

0.00186 

0.1*706 

P. 

1.37i^  X  106 

1.1*30  X  106 

1.377  X  106 

1.66  X  106 

where  Pi^  expresses  the  power  of  the  current  ratio  which  modiilates  P2,  the  zero 
time  value  of  X.  The  second  equality  was  derived  from  the  first  by  substitu¬ 
tion  of  Eq.  (1). 


The  characteristic  width  of  the  Gaussian  is  given  by 


^  =  P, 


P3t/T 


(H) 


where  P3  again  expresses  a  power  of  the  current  ratio,  this  time  modulating 
P]^,  the  initial  Gaussian  width.  (W  is  analogous  to  the  variance  of  the  normal 
distribution. ) 


P5  is  the  calculated,  zero-time  value  for  the  decaying  current  wave  form, 
Thou^  not  a  measured  quantity,  it  does  permit  use  of  an  analytic  formulation 
which  gives  current 

l(t)  =  P^  exp  (-t/x),  (5) 

accurately  representing  the  measured  current  when  the  armature  is  near  the 
probes.  <  '  ' 

P2  and  ?!  are  continuously  modulated  by  the  factors  defined.  This  repre¬ 
sentation  was  chosen  primarily  to  allow  compeurison  between  experiment  and 
models  whose  arc  length  varies  as  a  jxjwer  of  the  inverse  of  the  ratio,  instan¬ 
taneous  current  to  peak  current.  Note  that  Xp(t)  is  given  by  Eq.  (l).  The 
values  for  t  were  not  used  as  fitting  parameters,  but  were  taken  directly 
from  the  data.  In  this  way  the  time  dependence  of  the  current  was  forced  to 
match  the  experimental  conditions. 

After  numerous  iterations  of  the  five  parameters  listed  above,  minimiz¬ 
ing  the  differences  between  the  measured  and  calculated  B-dot  signals,  we 
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Parametric  Description  of  Armature  Current  Density  Function 


observed  the  parameters  shown  in  Table  IV,  All  of  the  parameters  except  Pi+ 
were  very  sensitive  to  the  fitting  procedure.  Since  Pp  is  not  a  large  number,  it 
is  evident  that  its  modulator,  P)^,  would  not  be  a  critical  parameter  in  deter¬ 
mining  the  calculated  B-dot  signal.  The  Gaussian  width,  P^,  and  its  modu¬ 
lator,  Pp,  were  sensitive  fitting  parameters  and  their  values  for  all  shots  are 
quite  similar. 

The  parameters  above,  together  with  the  positlon-versus-time  equation, 
completely  define  J  for  all  points  in  space  at  all  times.  The  results  can  be 
expected  to  describe  J  accurately  only  In  the  regions  of  time  when  some  signal 
is  present  on  the  B-dot  probes.  Beca\ise  J  is  a  function  of  two  variables,  we 
are  able  to  display  the  ciorrent  density  in  two  different  ways.  The  current 
density  J  may  be  plotted  either  as  a  fmiction  of  position  at  a  fixed  time  or 
as  a  function  of  time  at  a  fixed  position.  The  former  is  the  so-called  snap¬ 
shot  view  of  the  armatvire  current i  the  latter  is  representative  of  the  arma¬ 
ture  "flying  by"  a  point  on  the  barrel.  Figure  18  graphically  displays  the 
spatial  distribution  of  J  for  four  different  times  dxrring  Shot  96.  The  times 
have  been  arbitrarily  chosen  to  correspond  to  100  kA  Intervals  as  the  current 
decays  and  the  projectile  is  accelerated  down  the  bore.  It  is  evident  from 
these  curves  that  the  arc  length  is  between  20  and  ^0  cm.  Therefore,  the  peak 
current  density  does  not  exceed  150  MA  per  square  meter. 

b.  Arc  Length 

Since  the  arc  length  is  not  defined  for  a  Gaussian  function,  we  have 
arbitrarily  chosen  to  define  the  length  as  the  distance  from  the  rear  of  the 
projectile  to  a  point  in  space  at  which  the  c\irrent  density  has  fallen  to  one 
tenth  its  peak  value.  For  computational  purposes,  we  actually  Include  J  where 
it  was  larger  than  one  percent  of  the  maximum.  This  corresponds  to  a  point 
three  Gaussian  widths,  W(t),  behind  the  maximum  value.  Figure  19  displays  the 
"arc  length,"  calculated  in  this  manner,  versus  total  current  in  the  armature 
for  each  of  the  four  shots.  While  the  cutoff  point  of  ovir  distribution  is 
arbitrary,  the  trends  in  time  and  current  variation  are  real.  We  see  from 
this  plot  that  the  armature  length  scales  nearly  as  the  inverse  of  the  cur¬ 
rent. 


c.  Comparison  of  Current  Density  and  Light  Intensity 

Figures  20  and  21  show  plots  of  the  current  density,  J(t)  vers^ls  time 
for  the  fixed  locations  of  each  of  the  probes.  We  have  also  arbitrarily  scaled 
the  li^t  intensity  measxirements  to  the  peak  value  of  J  as  it  flies  by  each 
probe  location  and  plotted  it  for  comparison.  It  is  evident  that  luminous  ma¬ 
terial  not  only  outruns  the  projectile  in  Shot  96  but  also  lags  behind  the 
current  carrying  portion  of  the  armatvire  in  both  shots. 

Fine  structure  in  the  light  output  is  also  evident  in  both  shots,  particu¬ 
larly  at  later  times.  This  type  of  structure  was  previously  observed  in  data 
taken  from  the  BRL  railgun  and  can  possibly  be  explained  by  the  existence  of 
Raleigh-Taylor  instabilities  in  the  arc, 9 


9powell,  John  D. ,  "interchange  Instability  in  Railgun  Arcs,"  Phys.  Rev.  A, 
3^,  3262  (1986).  See  also  Ballistic  Research  Laboratory  Report  No.  BRL-TR- 
2TT6,  January  198'?’. 
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=  Shot  95 
=  Shot  96 
=  Shot  97 


The  length  was  defined  as  the  distance 
from  the  projectile  to  the  point  at  which 
J  reached  one- tenth  its  maximum  value. 
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Figure  21.  Comparison  of  Current  Density  and  Light  Emission;  Shot  96 
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III.  MODEL  CALCULATIONS 


In  the  earlier  parts  of  this  report  data  obtained  from  several  shots  on 
CHECMATE  have  been  discussed  in  some  detail.  It  is  obviously  desirable  to  ob¬ 
tain  in  addition  some  estimate  of  proi>erties  of  the  arc  which  csuinot  be  meas¬ 
ured  directly.  Such  properties  include,  for  exaniple,  the  arc  mass,  the  tem¬ 
perature  as  a  function  of  position  and  time,  and  the  extent  of  ionization  in 
the  arc.  In  this  section  of  the  report  we  discuss  some  calculations  that  we 
have  undertaken  to  determine  these  properties. 

The  model  used  for  the  calculations  is  basically  the  one-dtmenslonal, 
steady-state  model  previously  developed  by  Powell  and  Batteh.^  The  model  will 
be  applied  at  each  Instant  of  time  tmder  the  assumption  that  time-dependent 
changes  in  the  arc  properties  occur  in  a  quasl-statlc  manner.  This  assumption 
is  felt  to  be  reasonable  for  purposes  of  analyzing  the  data  discussed  here 
since  the  arc  is  rapidly  formed  and  esqperlmental  measurements  are  made  well 
after  the  acceleration  process  is  underway.  If,  Instead,  one  wished  to  study 
arc  formation  or,  say,  the  collision  of  the  arc  with  a  projectile,  then  highly 
time-dependent  effects  are  in5)ortant  and  must  be  accounted  for.  These  effects 
have  been  studied  by  us^®  and  others^^  in  previous  work.  There  are,  of  course, 
additional  effects  such  as  instabilities  and  turbulence,  not  yet  investigated 
in  great  detail,  and  these  effects  are  not  Included  in  the  model  or  calcula¬ 
tions. 


A  number  of  modifications  have  been  made  to  the  steady-state  model  both 
to  Inqprove  the  calculations  and  to  make  best  use  of  the  experimental  results. 
The  most  significant  change  is  that  we  have  omitted  entirely  the  theoretical 
calculation  of  temperature,  a  calculation  previously  based  on  the  assumption 
that  ohmic  heating  in  the  arc  is  baleuiced  at  any  time  by  radiative  cooling. 
Instead,  the  temperature  within  the  arc  is  inferred  from  the  measured  (fitted) 
current  density  as  a  function  of  position  and  the  measured  muzzle  voltage. 

This  modification  has  been  introduced  because  theoretical  temperature  calcula¬ 
tions,  particularly  in  one  dimension,  are  not  felt  to  be  reliable.  Not  only 
are  there  numerical  difficulties,  but  there  may  also  be  additional  energy  loss/ 
gain  mechanisms  not  accounted  for  and  difficult  to  quantify.  These  difficul¬ 
ties  are  avoided  by  simply  using  the  experimental  data  directly.  For  purposes 
of  comparison,  however,  we  have  vindertalen  scane  purely  theoretical  calculations 
in  which  temperature  is  calculated  via  energy- conservation  considerations, 
and  these  results  will  also  be  discussed  in  the  ensuing  pages. 

Two  other  minor  modifications  have  also  been  Introduced.  First,  the 
electrical  conductivity  in  the  model  has  been  extended  to  Include  not  only 
electron-ion  scattering  (Spitzer),  but  electron-neutral  scattering  as  well. 

The  latter  contribution  was  included  when  arc  temperatures  were  found  to  be 


^^owell,  John  D.  and  Jamison,  Keith  A.,  "One-Dimensional,  Time-Dependent 
Model  for  Rallgun  Arcs,"  Ballistic  Research  Laboratory  Report  No.  BRL-TR- 
2TT9,  Feb  87.  See  also  Jamison,  K.A.,  Burden,  H.S.,  Powell,  J.D.,  and 
Marquez-Reines,  M. ,  "Non-Steady  Phenomena  in  Railgun  Plasma  Armatures," 

IEEE  Trans.  Magn.  MAG-22,  15^6  (1986). 

^^Batteh,  Jad  H.  and  Rolader,  Glenn  E. ,  "Modeling  of  Transient  Effects  in  Rail- 
gun  Plasma  Armatures,"  Ballistic  Research  Laboratory  Contract  Report  No. 
BRL-CR-56T,  March  1987. 


31 


sufficiently  low  that  this  scattering  mechanism  might  in  some  cases  he  impor¬ 
tant.  The  second  change  in  the  previously  used  model  is  that  we  have  attempt¬ 
ed  to  account  very  approximately  for  the  finite  rail  height,  h  (see  Fig.  3). 

In  previous  work  the  rails  were  assumed  to  he  infinitely  high,  hut  this  assump¬ 
tion  is  known  to  overestimate  substantially  the  acceleration  and  pressure 
of  the  arc.  To  account  for  the  finite  rail  height  we  have  \ised  Batteh's 
approximation^^  in  which  the  magnetic  permeahillty  is  simply  replaced  hy  an 
effective  value. 


In  Sec.  III.l  we  outline  the  governing  equations  for  the  model,  discxjss- 
ing  in  particular  the  changes  noted  above.  In  Sec.  III. 2  we  present  results 
of  the  data  analysis,  indicating  properties  of  the  arc  not  directly  measured, 
and  compare  the  results  with  those  obtained  from  some  purely  theoretical  work. 
We  also  indicate  reasons  for  the  lack  of  agreement  and  suggest  ways  in  which 
future  models  can  he  Improved. 

1.  Governing  Equations 

We  take  directly  from  the  experimental  data  at  each  instant  of  time  values 
of  the  nuzzle  voltage  Vjjj,  the  arc  length  the  current  i,  the  projectile 
location  down  the  hore  Xp,  and  the  projectile  mass  mp.  We  also  essentially 
take  as  given  the  experimentally  determined  current-density  profile, 

1  , 

J  =  (JT  e  2W^  (6) 

N 

where  values  of  W  and  xq  are  also  provided  by  the  experiment  (see  Pig  IT). 

In  general,  however,  theoretical  expressions  for  the  arc  conductivity  vanish 
at  the  hack  of  the  arc  where  the  density  vanishes  so  that  it  is  convenient  to 
replace  Eq.  (6)  near  the  hack  of  the  arc  with  a  function  which  vanishes  at 
X  =  Xp  -  £a,.  For  simplicity  we  choose  a  quadratic  function  to  represent 
J  in  this  region.  In  the  notation  of  previous  model  calcixlations^  where  C 
is  a  dimensionless  distance  defined  hy 


C  = 


X  -  X  +  Jl 

_ E _ a 


we  therefore  have 


(T) 


®  »  for  o  (8) 

J(C)  =  C  +  Cp  e®  ,  for  C  <  C*  . 

Note  that  in  this  notation  C  =  0(x  =  Xp-£a,^  corresponds  to  the  trailing 
edge  of  the  arc  and  C  =  1  (x  =  X  )  corresponds  to  the  projectile  surface. 

The  value  of  C*  is  a 'hltrary,  hSt  should  he  fairly  small,  and  we  have  consist¬ 
ently  chosen  =  0.1.  The  constants  \Jq  and  Cq  in  Eq.  (8)  are  given  hy 

^^Batteh,  Jad  H.,  "Momentum  Equation  for  Arc-Driven  Rail  Guns,"  J.  Appl.  Phys. 
36,  3162  (196^).  See  also,  "Arc-Dynamic  Calculations  in  the  Rail  Gun," 
Ballistic  Research  Laboratory  Contract  Report  No.  ARBRL-CR-00521,  Nov  83. 


32 


and 


(9) 


«0  '  *"'a 


x„  -  X  +  £ 
r  =  O  p  a 
^0  £ 


the  constants  C^,  C2  and  C'  follow  from  the  requirement  that  J  and  be 
continuous  at  5*  and  from  the  current-conservation  relation 


/  Jd5  =  j/i 


(10) 


Here  ,1  is  the  rail  current  per  unit  height.  We  easily  determine  from  these 
conditions  that 


c- .  ^ (u) 

N  3j  Sj  j  '2^  ^  V,  ^ " 


2R. 


C,  = 


1  C-r» 


Z' 

N 


(12) 


and 


K 


2  C'r* 


N"  N 

In  these  expressions 


C' 


(13) 


2  w7  ® 


(lU) 


and  erf  (x)  is  the  error  function 


?  ^  -u® 

erf(x)  =  —  /  e  du  , 


(15) 


Once  the  current  density  has  been  obtained,  expressions  for  the  magnetic- 
induction  field  B  and  pressure  P  follow  immediately  from  results  derived  in 
Ref.  2.  In  particular,  we  have 


B( c)  =  uj  -  u£  J  j( , 


(16) 
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which  yields  after  employing  (8)  and  integrating 


B(C)  =  Uj  -  (C^CV2  +  C^CVS)  ,  for  C  <  C* 

and 


(17) 


B(C)  =  uj  -  Uilg^  (C^C*V2  +  C^K*^/3) 

-  terf  t-^) 

C-  ^2'  J2  W- 

N  0 


C 

erf  (  C  >  5* 


If  we  further  assume  that  the  projectile  mass  is  large  compared  to  the  arc 
mass,  as  we  will  find  to  be  the  case  in  all  our  analysis,  then  the  hydro¬ 
static  pressure  P  is  given  in  terms  of  the  B  field  by 


P 


(18) 


In  Eqs.  (l6)-(l8),  u  is  the  standard,  free-space,  magnetic  permeability. 
These  equations,  however,  have  been  derived  under  the  assumption  that  the 
rails  are  infinitely  high,  an  approximation  that  is  known  to  overestimate 
the  pressure  and  accelerating  force  of  the  arc.^  Batteh^^  has  eirgued  that 
approximate  account  can  be  taken  of  finite  rail  height  provided  u  in  these 
eqxiations  is  replaced  by  an  effective  value  u*»  In  a  tedious  calculation 
he  has  worked  out  the  appropriate  value  and  finds 


U*  "  fS  tan  ^  (S)  +  log  (l  +  1/S^)  - log  (l  +  S^)]  (19) 


where  S  =  h/w,  i.e.,  the  ratio  of  the  rail  height  to  rail  separation.  We 
will  employ  this  approximation  in  the  remainder  of  our  analysis. 


The  position-dependent  conductivity  within  the  arc  can  also  he  Inferred 
from  the  experimental  data,  namely,  the  current  density  in  (6)  and  the  measured 
muzzle  voltage.  It  has  been  notedi*^  that  in  the  steady  state  J/o  is  indepen¬ 
dent  of  position  in  the  arc  and  is  given  by 


J/o  =  Vp/w  .  (20) 

Here  o  is  the  arc  conductivity  and  Vp  is  the  potential  drop  across  the  plasma. 
Consequently,  we  have 

a  =  wJ/Vp  .  (21) 

3i< 


(22) 


In  terms  of  the  measured  voltage  Vjn,  Vp  Is  given  by 

V  =  V  -  V 
p  m  c 

where  V^,  is  the  total  contact  potential  at  the  two  rail-arc  interfaces.  Un¬ 
fortunately,  values  of  are  not  known,  so  we  have  analyzed  the  data  for 
various  reasonable  estimates. 

The  remaining  properties  of  the  arc  that  must  be  determined  are  the 
ten5)erature  T,  the  concentration  of  singly  and  doubly  charged  ions  and  X2# 
the  density  p,  and  finally  the  arc  mass  m^.  To  determine  the  first  three  of 
these  quantities  we  invoke  the  theoretical  expression  for  the  arc  conductivity 
along  with  the  Saha  equations  which  predict  the  ionization  state  as  a  function 
of  temperature  and  presstire.  For  the  conductivity  we  assume^^ 


1 


o 


(23) 


where  og  is  the  conductivity  of  a 
slightly  ionized  gas.  For  oh  the 


2.3  X  10~^  Yg 
Z  log  A 


t3/2 


completely  ionized  gas  and  ol  that  for  a 
standard  Spltzer  e;)^re8sion  is  used 


(2li) 


where  Z  is 


Z  = 


+ 

+  2x2 


(25) 


and  ye  constant  given  by  0.58T  for  a  singly  ionized  gas  and  by  0.683  for 
a  doubly  ionized  gas.  For  arbitrary  ionization  we  have  assiuned  for  sliipll- 
city  a  linear  relation 


Y„  =  0.096  Z  +  O.i+91 


(26) 


The  shielding  parameter  A  in  Eq.  (24)  can  be  written  in  terms  of  T,  P,  x^  and 
X2  as 

1.23  X  10^  (1  +  +  ^XjjkjT  *5 

* - Z - P(x^  .  2X5) - 

where  kg  is  Boltzmann's  constant. 


For  the  case  of  slight  ionization  we  use  the  relation 


^3cambel,  A.B.,  Plasma  Physics  and  Magnetofluldmechanics  (Mc-Graw-Hill,  New 
York,  1963),  Chap.  T. 
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(28) 


U.5  X  10“^^  (x^  +  2x2) 

'’l  ~  (1  -  x^  -  Xg)  Q  T^'S 

where  Q  is  the  electron-neutral  scattering  cross  section.  Values  of  Q  are 
generally  not  known  in  all  temperature  ranges  for  materials  of  Interest  so  we 
have  assumed  a  constant  value  of  10“ m*.  At  most  temperatures  a  is  not  a 
particularly  sensitive  function  of  Q. 

The  ion  concentrations  x^  and  Xg  follow  from  solutions  to  the  Saha  eq- 
viatlons  as  discussed  in  Ref.  2.  Thus,  we  have  the  relations 

x^  (x^  +  2Xp) 

Ti-v,  -  x^ni  -f  X, .  gx,)  -  ''1 

and  ( 29 ) 


Xg  +  2X2) 

x^(l  +  +  2x2) 


=  K2  (T,P)  . 


Calculation  of  the  functions  and  has  been  dlscvissed  elsewhere^  and  will 
not  be  discussed  here.  We  should  point  out,  however,  that  these  functions 
depend  only  on  the  temperature  and  presstare  and  on  the  the  properties  of  the 
material  which  constitutes  the  arc,  l.e.,  the  ionization  laotentials  and 
electronic  partition  functions. 

It  is  now  apparent  that  Eqs.  (23)  and  (29)  along  with  the  supplementary 
relations  discussed  and  the  experimentally  inferred  conductivity  can  be  solved 
simultaneously  to  yield  values  of  T,  x^^,  and  Xg  at  all  points  of  the  arc.  This 

solution  must  be  carried  out  numerically.  Once  the  values  are  known,  however, 
the  density  follows  from  the  ideal-gas  equation  of  state  assiosed  in  all 
o\ir  previous  work,  namely. 


P 


m^ 


(1  +  X, 


P 

+  2x2)kgT 


(30) 


where  mQ  is  the  atomic  mass  of  the  ions  or  neutrals  making  up  the  arc.  Fin¬ 
ally,  the  arc  mass  follows  directly  ftrom  the  calculated  density  and  measured 
arc  length  via  the  relation  - 

1 

m  =  h  w  J  p(c)dC  (31) 

a  a  ^ 


where  h  is  the  rail  height. 
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The  foregoing  analysis  is  sufficient  to  determine  all  the  arc  properties 
without  our  having  to  employ  any  relationship  relating  to  energy  conserva¬ 
tion.  We  now  tiirn  to  a  discussion  of  some  of  the  calculations. 

2.  Results  of  Calculations 

We  have  xmdertaken  numeroTis  calculations  using  the  formalism  of  the 
preceding  section  and  discuss  here  only  some  of  the  more  typical  results. 

The  calculations  fall  hasically  into  two  groups.  In  the  first,  certain  aver¬ 
age  or  position-independent  properties  of  the  arc  were  computed;  in  the  sec¬ 
ond,  spatially  dependent  quantities  were  calculated  at  a  given  point  in  time. 

All  data  were  taken  from  three  of  the  shots  on  CHECMATE  discussed  previously, 
those  numbered  95,  96  and  9T.  The  projectile  mass  for  Shots  95  and  97  was 
233g,  while  that  for  Shot  96  was  120g. 

The  material  which  constitutes  the  arc  was  not,  of  course,  known  in  any 
detail  and  most  of  o\ir  results  are  based  on  the  assumption  that  it  consists 
of  carbon  vapor.  Carbon  was  chosen  for  most  of  the  calculations  for  reasons 
which  will  be  explained  when  the  analysis  is  presented.  We  have,  however, 
undertaken  some  computations  for  copper  and  aluminiim-vapor  arcs  and  these 
results  are  briefly  discussed  primarily  to  indicate  some  of  the  more  salient 
differences.  We  have  also  assumed  in  most  calculations  that  no  contact  poten¬ 
tial  existed  at  the  rail-arc  Interfaces.  Thus,  the  total  measured  muzzle  voltage 
in  Eq,  (22)  is  assumed  to  be  a  bulk  voltage  drop  across  the  plasma.  We  do, 
however.  Indicate  in  some  of  the  later  discussion  what  effect  a  nonzero  contact 
potential  should  have  and  how  it  can  bring  experimental  and  theoretical  results 
into  better  agreement. 

Shown  in  Fig.  22  is  the  calculated  armature  mass  for  the  carbon  arc  as 
a  function  of  time  for  each  of  the  three  shots.  For  Shots  95  and  97  the  mass 
appears  to  remain  nearly  constant  at  about  1.5g  for  the  entire  acceleration 
time.  Thus,  if  a  significant  amount  of  ablation  is  occurring  during  the  shot, 
the  resulting  mass  does  not  appear  to  be  entrained  in  the  cu*c.  This  point 
will  be  discussed  further  when  a  comparison  with  some  pxarely  theoretical  re¬ 
sults  is  made.  For  Shot  96  the  arc  mass  is  sli^tly  smaller  than  for  the 
other  two  cases.  One  possible  explanation  for  the  slightly  smaller  mass  in 
Shot  96  is  that  the  injection  velocity  for  this  case  (smaller  projectile  mass) 
is  higher  than  that  for  Shots  95  and  97  and  the  dwell  time  at  any  given  point 
on  the  rails  is  also  smaller.  Consequently,  less  time  is  available  for  rail 
and  insulator  material  to  become  vaporized  and  entrained  in  the  arc  if  in 
fact  such  an  effect  occurs  at  very  early  times  before  the  nearly  steady  values, 
seen  in  Fig.  22,  are  achieved.  It  should  also  be  noted  that  there  are  fairly 
significant  drops  in  the  armature  mass  at  about  0,7  and  0.9  milliseconds. 

These  times  are  just  subsequent  to  the  times  at  which  light-output  data  from 
Shot  96  indicated  luminous  material  being  blown  pest  the  projectile  (see  Sec. 

II, 2, d).  At  each  time,  however,  the  equilibri\an  value  of  the  mass  appears 
to  be  quickly  restored. 

We  have  also  computed  as  a  function  of  the  pulsed  c\irrent  i  both  mean 
and  peak  temperatures  within  the  arc  as  well  as  the  mean  concentration  of 
single  ions  <x^>.  The  peak  temperature  occurs  just  to  the  left  of  the  pro¬ 
jectile  sirrface  as  will  be  seen  when  some  spatially  dependent  results  are 
presented.  These  results  are  shown  in  Figs.  23-25  and  exhibit  largely 
self-explanatory  behavior.  Both  temperature  profiles  increase  with  increasing 


1.60  r  “  =  SHOT  95 
0=  SHOT  96 


Calculated  Armature  Mass  Versus  Time  of  the  Carbon  Arc 


SHOT  95 


o 

O 

o 

O 

o 

o 

o 

O 

o 

O 

o 

o 

o 

O 

o 

o 

o 

o 

o 

00 

<£> 

CM 

o 

ro 

CNI 

<N 

CM 

CM 

CM 

3aniVd3dmi 


39 


Figure  23.  Peak  Temperature  Versus  Current  for  Carbon  Arc. 
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Figure  25.  Average  Ion  Concentration  Versus  Current  for  Carbon  Arc. 


current,  the  increase  resulting  primarily  from  enhanced  ohmic  heating  at 
higher  currents.  It  is  of  interest  to  note  that  very  approximate  scaling  re~ 
lations,3  derived  under  the  assumption  that  v  «  2^,  indicate  that  temperature 
shovild  increase  as  The  graphs  in  Figs.  23  and  2k  can  be  seen  to  obey 

roughly  such  a  relation. 

The  mean  concentration  of  single  ions  can  also  be  seen  to  increase  with 
increasing  current.  Basically,  this  behavior  results  from  con^jetition  between 
the  increasing  temperature,  which  enhances  ionization,  and  the  increasing 
pressure,  which  reduces  it.  The  ionization  degree,  however,  is  a  very  sensi¬ 
tive  function  of  temi^eratiare,  and  a  much  less  sensitive  function  of  pressure. 
The  concentration  of  doubly  charged  ions  was  also  computed  but,  at  these  rela¬ 
tively  low  temperatures,  their  values  were  too  small  to  be  evident  on  the 
graph. 


It  is  also  of  interest  to  ask  to  what  extent  the  results  of  the  data 
analysis  are  in  agreement  with  results  from  purely  theoretical  calculations. 

We  have  therefore  exercised  our  one-dimensional  model  to  undertake  a  number 
of  calculations  pertinent  to  the  experimental  results.  For  these  model  cal- 
c\ilations  it  is  necessary  to  specify  only  the  arc  length  SLg^,  the  projectile 
mass  ffip,  and  the  current  I.  All  remaining  qmntltles  are  then  calculated  as 
a  function  of  position  within  the  arc.  In  particular,  temperature  is  calcu¬ 
lated  by  the  energy-conservation  assvimption  that  ohmic  heating  is  balanced  at 
any  time  by  radiation.  Radiation  losses  are  represented  in  the  manner  dis¬ 
cussed  in  Ref.  2  except  that  here  we  have  scaled  the  resulting  temperatures 
down  by  the  factor  (Ag/A/^)^  where  Ag  is  the  surface  area  of  the  ends  of  the 
arc  and  A;^  is  the  surface  area  of  the  total  arc.  This  scaling  was  performed 
because  a  truly  one-dimensional  model,  which  accounts  for  only  losses  at  the 
ends,  is  known  to  overestimate  substantially  the  arc  temperature.  More  sop¬ 
histicated  approximations,  accoionting  for  losses  from  the  sides,  have  also 
been  imdertaken  but  this  simple  scaling  procedture  appears  to  be  nearly  as  ac¬ 
curate  as  any  of  them  and  is  very  simple  to  employ. 

Shown  in  Fig.  26,  Curve  A,  and  in  Table  V  are  the  theoretical  results^ 
for  the  temperature  as  a  function  of  position  within  the  arc  using  as  input 
the  parameters:  I  =  903  kA,  A  =29  cm,  and  m  =  0.233  kg.  These  parameters 
were  taken  directly  from  the  exjjerimental  data^for  Shot  95  at  the  time  t  = 

0.89  ns.  The  theoretical  temperatvire  profile  is  similar  to  that  we  have  ob¬ 
served  in  previous  model  calculations.  The  temperature  peak  occurs  Just  to 
the  left  of  the  projectile  and  steep  gradients  exist  near  the  projectile  sur¬ 
face.  The  steep  gradients  result  from  the  high  density  and  consequent  small 
radiation  mean  free  path  in  this  region.  In  essence,  photons  are  "trapped" 
behind  the  projectile  and  can  rpdlate  outward  only  when  they  are  within  a 
mean  free  path  or  so  of  the  boundary. 

In  the  theoretical  models  no  account  is  taken  of  any  contact  ixstentlal 
at  the  rail-arc  interfaces  so  that  only  the  potential  drop  across  the  plasma, 

Vp,  is  calculated.  This  potential  is  equivalent  to  the  experimentally  mea¬ 
sured  nuzzle  voltage  only  if  negligible  contact  potential  really  exists. 
Otherwise,  the  contact  potential  should  be  subtracted  from  the  muzzle  voltage 
^m  before  a  comparison  with  the  theoretical  results  is  xmdertaken.  Since 
numerical  values  of  the  contact  potential  are  not  known,  we  have  carried  out 
the  data  analysis  in  Fig.  26  for  three  different  values  of  V^,  namely,  0,  T5 
volts  and  l67  volts,  shown  in  ciirves  B,  C  and  D,  respectively.  The  last  value 


k2 


vas  chosen  so  that  the  potential  drop  across  the  plasma  Vp  was  the  same  as 
the  value  predicted  from  the  theory. 


TABLE  V.  Con^jarison  of  Theoretical  and  Data  Analysis  Results 


A 

B 

C 

D 

m 

a 

0.95  s 

1.U6  g 

1.25  g 

0.87  g 

V 

m 

- 

287V 

287V 

287V 

V 

c 

- 

0 

75V 

167V 

V 

P 

120V 

287V 

212V 

120V 

<T> 

2.77  X  10^ 

2.03  X  10^ 

2.19  X  10^ 

2.72  X 

H 

0 

T 

MAX 

3.71*  X  10^ 

2.52  X  10^ 

2.76  X  10^ 

3.62  X 

lo'* 

<X^> 

0,lh 

0.U8 

0.59 

0.77 

<Xp> 

0,03h 

0.003 

0,00k 

0.025 

It  Is  evident  from  Fig.  26,  curves  A  and  B,  that  there  is  substantial 
disparity  between  theory  and  experiment  if  one  assiimes  that  no  contact  poten- 
tial  exists  at  the  rail-arc  interfaces;  at  some  points  the  two  temperatures 
differ  by  as  much  as  50%.  As  V^.  increases,  however,  as  shown  in  curves  C  and 
D,  the  experimental  temperatiires  also  increase  and  approach  the  theoretical 
profile.  The  increase  in  T  with  increasing  Vj.  results  from  the  consequent 
smaller  values  of  Vp.  These  smaller  values  can  result  only  if  the  ''onducti- 
vity,  and  thus  the  tenqjerature ,  also  increases.  It  seems  probable  u.:at  the 
contact  potential  lies  somewhere  between  the  extreme  values  of  0  and  l67  volts 
so  that  a  realistic  curve  from  the  data  analysis,  accounting  accurately  for 
the  contact  potential,  would  still  lie  somewhat  below  the  theoretical  profile. 
If  so,  it  appears  that  the  theory  probably  overestimates  the  temperature  of 
the  arc  by  a  few  percent.  One  likely  source  of  error  in  the  model  is  the 
neglect  of  some  energy-loss  mechanism.  Recently,  it  has  been  proposed  that 
neutral  atoms  may  diffuse  out  the  back  of  the  arc,  while  mass  is  added  via 
eblatlon  in  such  a  way  that  the  toteJ.  arc  mass  renalns  nearly  constant.  If 
so,  the  neutrals  will  carry  energy  with  them  and  this  loss  mechanism  is  not 
presently  accounted  for  in  the  theory.  That  there  is  high  concentration  of 
neutrals  near  the  back  of  the  arc  can  be  seen  in  Fig.  27  in  which  the  neutral 
concentration,  averaged  over  the  back  25%  of  the  arc,  is  plotted  as  a  func¬ 
tion  of  current.  The  data  used  to  ccmqnite  this  curve  are  identical  to  those 
used  in  the  computation  of  results  in  Figs.  22-2h.  It  seems  unlikely  that 
such  large  numbers  of  neutrals,  particularly  at  low  currents,  could  be  car¬ 
ried  along  with  the  arc  without  some  significant  loss  out  the  back. 


UU 


006*0 


00*^ 


♦ 

• 

o  ♦ 
o  » 
o  ♦ 
o  ♦ 

/" 

.*”»!> 

^  o 
^  *0  0 
♦  n®  » 

D«  “ 


♦  O  « 

♦  o  o 


10  to  rv. 
Q>  Q>  o 

I—  »—  ^- 

o  o  o 

XXX 
CO  to  (/) 

11  II  II 


o  Z 
LJ 
(T 
tx 

CO  3 

^  O 


in 

o 

m 

o 

m 

o 

m 

in 

CM 

o 

in 

CM 

00 

00 

00 

00 

• 

d 

• 

o 

d 

• 

o 

d 

d 

o 

N0llVyiN33N0D  IVdinBN 


45 


Figure  27,  Neutral  Concentration  Averaged  over  Trailing  Edge  of  Carbon  Arc 


Unfortunately,  we  cannot  determine  \mequi vocally  the  cause  for  the  dis¬ 
crepancy  between  the  theoretical  and  experimental  temperature  profiles.  No 
doubt  both  effects  discussed  in  the  preceding  paragraph  contribute  to  some 
extent  to  the  lack  of  agreement.  Crude  estimates  which  we  have  made  suggest 
that  the  energy  lost  via  neutral  diffusion  is  probably  fairly  small  (at  most 
of  the  order  of  10%  of  the  arc's  internal  energy)  but  these  studies  require 
more  analysis.  We  are  not  aware  of  any  reliable  estimates,  theoretical  or 
experimental,  of  the  contact  potential. 

We  should  also  point  out  that  it  wo\ild  obviously  be  more  desirable  to  com¬ 
pare  results  of  the  data  analysis  with  theoretical  calculations  obtained  from 
a  model  in  which  the  contact  potential  was  specifically  inclioded.  One  covild 
then  determine  to  what  extent  the  existence  of  such  a  potential  affected  the 
internal  properties  of  the  plasma.  Unfortunately,  no  such  model  exists.  It 
does  seem  reasonable,  however,  to  exi)ect  that  such  effects  would  not  be  too 
significant  since  the  potential  is  localized  near  the  rail  stirfaces.  Further, 
studies  are  nonetheless  appropriate  when  the  nature  and  magnitude  of  the  con¬ 
tact  potential  are  better  understood. 

We  have  also  carried  out  some  calculations  under  the  assumption  that  the 
arc  was  composed  of  aluminum  and  copper  ions.  A  summary  of  the  results  of 
these  calculations  is  Indicated  in  Figs.  28  and  29  in  which  we  have  plotted  as 
a  function  of  current  the  arc  mass  i%  and  the  mean  temperature  T  for  the  three 
different  types  of  arc  materials.  All  data  were  taken  from  Shot  95  so  that 
the  ctu*ves  for  carbon  were  shown  in  Pigs.  22  and  2U.  For  reference  the  first 
and  second  ionization  potentials  and  atomic  masses  of  the  materials  involved 
are  shown  in  Table  VI.  As  can  be  seen  in  Fig.  28,  the  masses  of  the  different 
arcs  are  substantially  different  and  scale  very  rotighly  with  the  atomic  mass. 
Furthermore,  it  is  noteworthy  that  the  mass  computed  for  the  aluminum  arc  kg) 
is  significantly  higher  than  was  the  mass  of  the  fuse  initially  blown  to  create 
the  arc  ( »  0.3g).  Since  no  other  source  of  aluminiim  was  present,  it  then  ap¬ 
pears  that  only  a  very  small  percentage  of  the  arc  actually  consists  of  alumi¬ 
num  vapor.  Similar  statements  cannot  be  made  about  the  other  materials.  How¬ 
ever,  since  copper  has  a  thermal  conductivity  higher  than  the  insulating  mater¬ 
ial,  and  is  therefore  harder  to  vaporize,  it  appears  likely  that  the  arc  is 
comixjsed  primarily  of  elements  from  the  insulators,  i.e.,  carbon,  hydrogen  and 
oxygen.  The  temperatures  of  the  various  arcs  are  not  nearly  so  different  as 
their  masses  but  in  general  increase  with  increasing  ionization  potential. 

The  result  can  be  vinderstood  if  it  is  remembered  that  all  three  arcs  have  a 
common  conductivity  and  pressure  and  thus  hi^er  temperatTires  are  needed  to 
ionize  the  materials  which  have  higher  ionization  potentials.  However,  the 
extent  of  ionization  is  a  sensitive  function  of  the  ten5)erature  so  fairly  small 
changes  in  this  parameter  mean  large  changes  In  ion  concentration. 


TABLE  VI,  Ionization  Potentials  and  Atomic  Masses  of  Arc  Materials 

I^(ev)  Igtev)  m^Camu) 


Material 


2000 


3anivy3d^N3i 
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Figure  29.  Mean  Arc  Temperature  Versus  Time  for  Different  Arc  Materials 


IV.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  diagnostic  measurements  and  analysis  imdertaken  here  have  led  to  a 
number  of  conclusions  concerning  plasma  arcs  in  large-bore  railguns.  These 
conclusions  may  be  summarized  as  follows.  First,  It  appears  that  the  total 
mass  of  the  arc  remains  nearly  constant  over  the  acceleration  time.  Such  an 
effect  could  occur  if  there  were  no  ablation  from  the  rails  or  insulators  and, 
in  addition,  no  mechanism  for  losing  mass  from  the  arc.  It  is  more  likely, 
however,  that  some  ablation  does  occxir,  particularly  from  the  insulating  ma¬ 
terials,  and  the  resulting  Increase  in  mass  is  con^nsated  for  by  some  loss 
mechanism.  Two  possible  mechanisms  are  diffusion  of  neutrals  out  the  back  of 
the  axe  and  losses  at  boundary  layers  whose  velocity  must  vanish  at  the  rail 
surface.^^  Both  these  problems  merit  further  study  as  do  efforts  to  quantify 
the  extent  of  ablation  as  a  function  of  the  other  arc  properties. 

Second,  it  appears  that  theoretically  derived  teii5>eratvires  are  still  sig¬ 
nificantly  higher  than  those  inferred  from  experiment.  This  observation  per¬ 
sists  even  though  the  one-dimensional,  theoretical  model  has  been  revised  to 
accomit  approximately  for  radiation  losses  from  the  sides  of  the  arc.  It  is 
possible  that  the  theoretical  results  fail  to  acco\ant  for  some  energy-loss 
mechanism  (such  as  would  acconqjany  the  mass  loss  described  above)  and  thus 
overestimate  the  temperature.  It  is  also  possible,  however,  that  there  is  a 
significant  contact  potential  at  the  rail-arc  interfaces  so  the  measured  muz¬ 
zle  voltages  are  much  higher  than  the  potential  drop  across  the  plasma.  As 
has  been  seen,  postulating  a  contact  potential  and  subtracting  its  value  from 
the  muzzle  voltage  before  undertaking  the  data  analysis  leads  to  substantially 
better  agreement  of  theory  and  experiment.  Undoubtedly,  both  these  factors 
are  present  and  the  extent  of  each  can  be  determined  only  by  further  study. 

Of  particular  Importeince  are  escperlmental  efforts  to  measure  the  contact 
potential  and  theoretical  studies  of  the  mass  diffusion  postulate. 

Finally,  it  does  not  appear  possible  to  gain  much  information  concerning 
the  composition  of  the  arc  from  the  work  underteJeen  here.  It  is  fairly  evi¬ 
dent  that  there  is  very  little  aluminiim  present  but  the  dominant  material  can¬ 
not  be  inferred.  Presumably,  the  arc  is  formed  at  very  early  times,  before 
its  apparently  steady  mass  is  achieved,  by  material  ablated  from  the  projec¬ 
tile,  rails,  or  insulators.  Experimental  efforts  to  determine  the  composition 
of  the  arc  would  obviously  be  worthwhile  and,  then,  theoretical  models  could 
be  revised  to  account  for  the  appropriate  mixture.  It  does  appear,  however, 
that  many  arc  properties  are  not  significantly  affected  as  the  arc  composi¬ 
tion  is  varied  among  several  reasonable  alternatives. 


Barber , 


John  P. , 


private  communication. 
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